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ABSTRACT 


A recently developed analytic result in acoustics, “Formulation IB,” is used to compute broad- 
band trailing edge noise from an unsteady surface pressure distribution on a thin airfoil in the time 
domain. This formulation is a new solution of the Ffowcs Williams-Hawkings equation with the 
loading source term, and has been shown in previous research to provide time domain predictions 
of broadband noise that are in excellent agreement with experimental results. Furthermore, this 
formulation lends itself readily to rotating reference frames and statistical analysis of broadband 
trailing edge noise. Formulation IB is used to calculate the far field noise radiated from the trailing 
edge of a NACA 0012 airfoil in low Mach number flows, by using both analytical and experimental 
data on the airfoil surface. The acoustic predictions are compared with analytical results and ex- 
perimental measurements that are available in the literature. Good agreement between predictions 
and measurements is obtained. 
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1. INTRODUCTION 


Trailing edge (TE) noise has been the subject of extensive research within the aeroacoustic 
community for decades, both experimentally and analytically. Areas of current research include 
the prediction of TE noise from rotating machinery and airframes. Research in the area of TE noise 
prediction has, in large part, been motivated by the desire to incorporate the results of TE noise 
analysis into a design methodology. The present work is similarly motivated, and the resulting 
formulation should lend itself well to an engineering design tool suite when aeroacoustics plays a 
role in the design. 

The literature abounds with various theoretical approaches to TE noise prediction. Howe [1] 
categorized the various theories of TE noise into three groups: 

i. Theories based on the Lighthill [2] acoustic analogy, e.g., Ffowcs Williams and Hall [3]. 

ii. Theories based on the solution of special problems approximated by the linearized hy- 
drodynamics equations, e.g., Amiet [4, 5] and Goldstein [6]. 

iii. Ad hoc models involving postulated source distributions for which strengths and types 
are empirically determined. 

The present work falls into the first category. “Formulation IB,” originally proposed in [7], is a 
solution of the loading source term of the Ffowcs Williams-Hawkings equation [8]. Such time domain 
methods allow a total decoupling of the acoustic signal from the aerodynamics. As such, these 
methods readily avail themselves to acoustic predictions that are based on input from experimental 
measurements or computational fluid dynamics (CFD) solutions. For example, Singer et al. [9], used 
a solution of the Ffowcs Williams-Hawkings equation to predict TE noise from sources that were 
modeled with CFD simulations. The acoustic formulation in their work is known as “Formulation 
1A” [10]. What distinguishes Formulation IB from this prior formulation is its relative simplicity, 
which makes it highly suitable for rotational reference frames and statistical analysis of broadband 
TE noise. 
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In [7], Formulation IB was applied to the prediction of far field noise due to incident turbulence 
on a NACA 0012 airfoil at tunnel speeds ranging from 40 m/s to 165 m/s, and the predictions 
were compared with the experimental results of Paterson and Amiet [11]. The time-dependent 
surface pressure required as input to Formulation IB was generated by stochastic modeling of the 
incident turbulence and approximation of the airfoil response with a result from thin airfoil theory 
[12]. Formulation IB was then used to predict the acoustic pressure as a function in time at a 
prescribed microphone location. The time domain results were then Fourier analyzed to determine 
the narrowband spectra of the far field noise. The far field spectra were found to be in excellent 
agreement with the frequency domain predictions and experimental measurements of Paterson and 
Amiet [11]. 

In the present work, the time domain approach described in [7] is used to predict far field 
radiation from the trailing edge of an airfoil. In the following section, Formulation IB is briefly 
reviewed for the case of a flat surface in a general nonunifornr motion. Some advantages of this new 
formulation relative to other solutions of the Ffowcs Williams-Hawkings equation are described. 

In Section 3, a model problem is considered in which a single- frequency surface pressure induces 
an acoustic source at the trailing edge of a flat plate in uniform motion. The unsteady surface 
pressure is an analytical result from thin airfoil theory that is taken from the work of Amiet [4, 
5, 12]. Two simple test cases are presented for validation purposes. The directivity of the tone 
induced by this surface pressure is examined for qualitative correctness. The results of a velocity 
scaling exercise are shown to be consistent with the results of Ffowcs Williams and Hall [3]. 

In Section 4, the surface pressure formulation introduced in Section 3 is extended by linear 
superposition to provide an analytic source model for broadband TE noise. This stochastically 
modeled surface pressure is used as input to Formulation IB to predict broadband TE noise from 
a NACA 0012 airfoil in a low-turbulence uniform mean flow. The surface pressure correlations 
that are required in the aerodynamic model are taken from two sources: an empirical flat plate 
formulation [13, 14] and experimental data [15]. The resulting calculations are compared with the 
acoustic predictions of Schlinker and Amiet [16] and the experimental measurements of Brooks and 
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Hodgson [17]. 


2. ACOUSTIC FORMULATION 

In [7], Formulation IB was derived for the case of a flat plate in a general planar motion. Only 
the result is presented here; in order to do so, the following discussion is required. 

Consider a flat, finite surface moving in the plane X 3 = 0 along a velocity vector V. The velocity 
vector V and the surface geometry are related to the coordinate axes as pictured in Fig. 1. Let 
/(xq, x '2 , t) denote a geometric function where / = 0 on the surface edge and / > 0 on the interior 
of the surface. Let v — V/ denote the unit normal that lies in the plane of the surface, is normal 
to the edge, and is directed inward (Fig. 1). Note that V need not be constant in space or time. 
The only stipulation on the velocity is that the motion of the surface is in the same plane as the 
surface. 

Let x = [xq, X 2 , X 3 U denote the position of an observer, and y = [yi,y2,0] T denote the position 
of a source point on the surface (Fig. 1). The spatial frame of reference is fixed to the undisturbed 
medium. The unsteady surface pressure p(y, r) gives rise to sound that initiates at source time r 
and radiates to the observer along r — x — y. This sound is denoted by p'{x. t). the perturbation 
pressure that arrives at the point (xi, £ 2 , 2 : 3 ) a t time t. Formulation IB predicts this sound by: 


4n p'(x,t) = 


f> 0 


( dp /dr — V dp/ ds ) cos 0 


eg r ( 1 — M r ) 


dS 

ret 


( 1 ) 


+ 


f 

p cos 6 

dS- 

ret *> 

f 

M v p cos 6 

h> 0 

_r 2 {l-M r )_ 

7=0 

r ( 1 — M r ) 


di , 

ret 


where M v and M r are the Mach numbers in the directions of i> and r, respectively. The quantity 
dp/ds is the pressure gradient in the direction of V, and V is the local magnitude of V. The 
subscript “ret” denotes evaluation at retarded time t — v/cq. 

The first and third integrals in Eq. (1) are far field terms, and the second integral is the near 
field term. As previously discussed in [7], the relative contributions of the terms in Eq. (1) are 
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considered, under the conditions of low Mach number and an observer in the acoustic far field, i.e., 


M<1, r»A (2) 

where A is a typical acoustic wavelength of interest. With respect to M and r, the surface far 
field integral, i.e., the first integral in Eq. (1), is proportional to 1/r, whereas the second and third 
integrals are proportional to 1/r 2 and M/r, respectively. Therefore, the far field surface integral 
dominates the signal under the conditions in Eq. (2). 

Perhaps most significant in regard to the form of Eq. (1) is that it is valid, as is, for rotating 
surfaces. Its predecessor, Formulation 1A [10], is significantly more complicated in its rotational 
form and cannot be approximated by only one surface integral in the far field. Such a significant 
simplification for far field calculations makes Formulation IB more suitable for statistical analysis 
of broadband sources for rotating surfaces. A statistical formulation based on Eq. (1) was derived 
in [7]. The focus of the current work is the application of Formulation IB to TE noise in the form 
of Eq. (1), as will be demonstrated in the following two sections. 


3. MODEL PROBLEM TRAILING EDGE TONE 

Any noise prediction made with Eq. (1) will be only as good as the input surface pressure p(y, r). 
The current thinking is that such time-dependent pressure data would result from experimental 
measurement or a CFD calculation. However, in this section, a simplified analytic expression is 
used for p(y, t ) to serve as a model problem. A result from thin airfoil theory [4] will be used 
to describe the unsteady surface pressure that is produced by the passage of a single frequency 
disturbance past the trailing edge of a slender airfoil. This simple surface pressure formulation will 
be extended to a broadband source model in the following section. 

3.1. SURFACE PRESSURE FROM THIN AIRFOIL THEORY 

The airfoil for this model problem is a rectangular flat plate in the plane x% = 0, undergoing 
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a uniform rectilinear motion, as in Fig. 2. The velocity vector V = [— U, 0, 0] T , where U is a 
constant subsonic speed. The plate’s surface and its boundary, / > 0, are defined by the rectangle 
{ — C < x\ <0} x {—b<x 2 < &}, with the trailing edge at x\ = 0. An unsteady pressure 
distribution is assumed on this surface, and is analytically prescribed from thin airfoil theory, as 
discussed below. 

Amiet [4] has proposed a formulation to model the response of an airfoil to the passage of a 
pressure disturbance over its trailing edge. This formulation, formally derived in [12], is based on 
the theory of a thin airfoil of infinite span and models the moving disturbance as stationary in the 
variable x\ — U c t , where U c is the convection speed of the disturbance. The induced pressure jump 
on the airfoil surface can be written 

AP( Xl ,t) =2P 0 g(x ll k c )e- ik ^- u ^ (3a) 

where k c = u/U c is the streamwise convective wave number, and Po is the amplitude of the distur- 
bance. The factor of two in Eq. (3a) indicates that the pressure is assumed to be antisymmetric 
between the upper and lower surfaces, and this expression thereby accounts for the pressure on 
both sides of the airfoil, i.e., the pressure jump. Note that Eq. (3a) differs from the general form 
for the pressure jump in [7] because the explicit term e tkcX ' in Eq. (3a) was incorporated into the 
transfer function g in [7]. The formulation in Eq. (3a) is used here for consistency with the TE 
noise research of Schlinker and Amiet [16]. 

The transfer function g(x i, k c ) is 

g{xi,k c ) = -1 + (1 + i)E*[-xi(k c + g(l + M))] (3b) 

where g = Mw / (3 2 U , f3 = y/1 — M 2 , and the function E* is given by 

E*(0 = f T~~ Y du = m-iS{0 (3c) 

Jo ( 2j7TU ) 2 

The quantities C(() and <S(£) are the Fresnel cosine and sine integrals, and will be evaluated 
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numerically by the formulas in [20]. The final representation for the unsteady surface pressure 
p(yi,r), assumed to be a real quantity, is 

p(yi,T) = S{-AP(yi,r)} (3d) 

The pressure jump is negative in Eq. (3d) because of the orientation of upper and lower surface 
normals in the derivation of the acoustic formulation [7] . 

Note that the transfer function in Eq. (3b) represents the effect of the induced surface pressure 
only, and neglects the effect of the incident pressure. The neglect of the incident pressure field effect 
is not of concern here, as this model problem is presented for illustrative purposes only. After the 
initial derivation of this induced pressure formulation [4, 12], Amiet later altered the formulation to 
include the effect of the incident pressure field [5] . The effects of both induced and incident surface 
pressure will be employed in the broadband formulation in Section 4. 

3.2. DIRECTIVITY CALCULATION 

Using Eqs. (3a)-(3d) as the input surface pressure in Eq. (1), the directivity of a single frequency 
source is now examined. The flat plate has a chord length C = 0.5 m, and a span 2 6 = 2.0 m. 
The flow speed U is determined by a freestream Mach number M = 0.2, with cq = 343 m/s. The 
disturbance amplitude Pq is taken as one percent of the dynamic head Qq = pqU 2 / 2. with p 0 = 1.23 
kg/m 3 , and the convection speed is taken to be U c = 0.8 U . The initial surface pressure p(y i,0), 
— C < y\ < 0. is shown in Fig. 3. This pressure profile represents the surface pressure over the 
entire span at observer time t = 0. Note, again, that the formulation in Eqs. (3a)-(3d) represents 
the induced surface pressure only. 

The radiated noise p'(x,t) is calculated at 360 equally spaced locations on a circular arc in the 
plane X 2 = 0. The radius of this arc emanates from the mid-span location on the trailing edge, as 
shown in Fig. 4. The arc trajectory (r, ip) is determined by r = 2 nr and 0 < ip < 2tt. The surface 
discretization is a uniform grid of 100 x 400 surface elements. The directivity is determined by the 
peak pressure amplitude calculated at each position, during one period T = 1/ f for a frequency of 
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2.5 kHz, with 128 timesteps in a period. Fig. 5 shows the results, in polar form, where the notation 
| \p' 1 1 is interpreted as 

lb'll = 0 ™^b'OM) I 

The upstream directivity of the major lobes is consistent with the research of previous authors, 
e.g., Singer, et al [8]. 

3.3. VELOCITY SCALING LAW 

Attention is now turned to the way in which the intensity of the far field noise, as predicted by 
Eq. (1), will scale as a function of velocity, when the surface pressure is described by Eqs. (3a)-(3d). 
A scaling law will be determined under the assumption that the acoustic source is noncompact, 
i.e., A C C. Furthermore, the observer is assumed to be in the acoustic and geometric far field, 
i.e., r> A and r >> C. respectively. 

Because scaling laws are typically determined for low Mach number flows [3, 19], the Mach 
number range of interest is 0.01 < M < 0.2. The surface pressure amplitude Po is one percent of 
the dynamic head. The physical dimensions of the plate are the same as in the above directivity 
problem. The observer is chosen at a distance of 10 m, directly above the trailing edge, i.e., 
x = [0,0, 10] T in meters. The calculations are performed on a 100 x 400 uniform surface grid. 

The surface pressure in Eqs. (3a)-(3d), with a frequency of 2.5 kHz, is used as input to Eq. (1) 
to predict the far field sound p'(x,t) to the observer. A separate calculation is run for each of 50 
equally spaced Mach numbers between 0.01 and 0.2. Each calculation is performed for one period 
with 128 timesteps. The average intensity I(x) of the acoustic signal at the observer x. assuming 
spherical spreading, is then calculated by 

= - [ T 

T Jo po co 

The average acoustic intensities for this test case, as a function of Mach number, are represented 
as circles in Fig. 6. The slope of these results on a log-log plot can be visually determined by 
observing their proximity to the dotted line with a slope of exactly five. This U 5 proportionality is 
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consistent with the result of Ffowcs Williams and Hall [3], as expected from the idealized conditions 
placed upon the calculations. 

4. BROADBAND TRAILING EDGE NOISE 

The analytic surface pressure in the previous section is extended to model a broadband trailing 
edge source on a slender airfoil at zero angle of attack. Following the approach of Schlinker and 
Amiet [16], the surface pressure correlations required as input are evaluated by flat plate theory and 
by experimental measurements. This broadband surface pressure is used as input to Formulation 
IB to predict far field radiation in the time domain. The results are Fourier analyzed and compared 
with experimental TE noise spectra [17]. 

4.1. EXPERIMENT DESCRIPTION 

The experiment modeled in this section is reported by Brooks and Hodgson [17]. A NACA 
0012 airfoil is placed between two plates at zero angle of attack in the test section of an open jet 
wind tunnel. A schematic of this experimental setup is shown in Fig. 7. Noise propagates from the 
test section into an anechoic chamber that is instrumented with microphones. 

The airfoil has a chord length of 0.6096 m and a span of 0.46 nr. The tunnel speeds of interest 
here are 38.6 m/s and 69.5 m/s. The chord-based Reynolds numbers are 1.57 million and 2.82 
million, respectively. Boundary layer tripping was applied at 15 percent chord downstream of 
the leading edge to ensure a spanwise uniform transition location and a fully developed turbulent 
boundary layer at the trailing edge. 

For radiated noise measurements, eight microphones are located in the plane perpendicular to 
the airfoil midspan. The presence of extraneous noise sources precludes direct measurement of TE 
noise by a single microphone. Therefore, to evaluate the TE noise, a cross-spectral analysis of pairs 
of microphones was employed in a manner consistent with the coherent output power method [20, 
21]. The microphone pictured in Fig. 7 represents the location for which the current predictions 
are made, at a distance of 1.22 m directly above the airfoil trailing edge. Note that a shear layer 


10 of 38 



forms downstream of the nozzle lip, between the airfoil and the microphone. Although both the 
directivity and the amplitude of the TE noise are affected by refraction through this shear layer, 
the corrections for the microphone at this location are small enough to ignore (see [17]). 

4.2. BROADBAND ANALYSIS 

For prediction purposes, the airfoil is modeled as a flat plate in order to evaluate the unsteady 
surface pressure with a broadband extension of the analytic formulation in Section 3. The airfoil 
geometry is oriented with respect to the coordinate axes as in Fig. 2, with { — C <aq<0}x{ — b < 
x 2 < b}, where C = 0.6096 m and 2b = 0.46 m. 

The surface pressure arises from boundary layer turbulence that is assumed to convect in a 
frozen pattern along the airfoil surface toward the trailing edge. Unlike the single frequency source 
in Section 3, the broadband nature of the surface pressure in the present case requires consideration 
of both chordwise and spanwise wave numbers, k\ and k 2 , respectively. Each Fourier component 
of this broadband surface pressure jump is associated with a wave number pair (ki, k 2 ) and can be 
written 

AP(h,k 2 -, Xl ,x 2 ,t) = 2P(A: 1 ,A: 2 ) ff (xi,A; 1 ,A: 2 )e-*[ fc 4^t)+fc 2 x3] (4) 

where P(ki 1 k 2 ) is the amplitude of the pressure jump. Because of the assumed convective nature of 
the turbulence, the notation for the chordwise wave number k\ will be replaced by k c to emphasize 
its dependence on the convection speed U c . In an exact sense, an infinite number of combinations 
of frequency and convection speed ratio u/U c would yield a given value of k c = lo/U c . However, the 
acoustically relevant structures in the turbulent boundary layer are assumed frozen with respect to 
a single convection speed that is taken as U c = 0.8 U. 

The complete broadband surface pressure jump AP(xi, X 2 , t) is obtained by summing all Fourier 
components in Eq. (4): 

/ oo roo 

/ f(Ae,* 2 )s(iiA,fe)e- i[ ‘ t( * , - p<, ) + ‘ ,I, ]<JMfc 2 (5) 

-oqJ — 00 

A straightforward approach for predicting the desired broadband far field measurements is to 
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use the real part of Eq. (5) as input to Formulation IB. This approach requires knowledge of a 
two-component surface pressure spectrum and a dual wave number transfer function g(xi, k Cl A^). 
However, because one of the objectives of the current work is to reproduce the results of Schlinker 
and Amiet [16] from a time domain perspective, an approach similar to theirs taken will be used 
to model the surface pressure. 

The analysis for the frequency-domain formulation in [16] comes from previous work [22] in 
which Amiet argues that only one spanwise wave number contributes to the sound detected by an 
observer in a given location. In particular, Amiet focuses on an observer in a spanwise symmetric 
location, for which only the zero spanwise wave number needs to be considered. An analogous time- 
domain result was derived in [7] . Although the analyses in [22] and [7] were presented to derive an 
acoustic formulation for incident turbulence noise, the results pertaining to spanwise wave numbers 
are sufficiently general to apply to the present trailing edge problem, with the only significant 
difference being the transfer function g{x\,k\,k 2 )- Therefore, applying the time-domain analysis 
in [7] to the present problem yields the following simplified expression for the surface pressure: 


AP{xi,t) =2ir P(k c ,0)g{xi,k c ,0)e tkc{ - xi Uci) dk c 

J — OO 

( 6 ) 

p{y,T) = 3?{-AP(t/i,t)} 


Note that Eq. (6) is valid as an acoustic source pressure for an observer in the plane X 2 = 0, 
and the k -2 integration in Eq. (5) is neglected. Furthermore, all three integrals in Eq. (1) will be 
evaluated for the predictions that follow, although the first integral is expected to dominate the 
signal. 

The unsteady surface pressure is numerically evaluated by an approach similar to that in [7]. 
Eq. (6) can be approximated by a truncated series in which the coefficients are chosen such that 
the autocorrelation and power spectrum of the series form a Fourier transform pair. This approach 
is consistent with stochastic modeling theory, e.g., [23]. The unsteady surface pressure jump in Eq. 
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(6) is then approximated by 


N 

AP(x u t) « 2-k ^2 A n e t(pn g{x 1 ,k c , n ,0)e-' lkc ’ n ( Xl - Uct ' > (7a) 

n=—N 

k c ,n = n Ak c , n = 0, ±1, ±2, . . . , ±N 
Ak c = k CjN /N 

The discrete surface pressure amplitudes {A n } are evaluated by 

A n = [$pp(A: Cj „,0) Afc c ]s (7b) 

where Qpp(k Cl ^ 2 ) is the two-component power spectral density (PSD) of the surface pressure. 
Schlinker and Amiet [16] argue that the required single wave number spectrum Qpp(k c ,0) can be 
evaluated by 

$pp(K,0) = hiw) S qq {u,0) (7c) 

where i 2 (cu) is the spanwise correlation length and S qq (u), Aa^) is the spanwise surface pressure 
correlation function. 

The phase angles are independent random variables uniformly distributed on [0, 2tt]. 

The transfer function in Eq. (3b) can be used for g(yi, k c ^ n , 0) with the following modification. As 
previously noted, Eqs. (3a) and (3b) represent the induced pressure jump. Amiet [5] has suggested 
that the incident pressure, i.e. , the pressure that results from turbulent eddies contacting the 
trailing edge, can be accounted for by the addition of an exponential convergence factor of the form 
e e kc x 1 , where e is a positive parameter. For — C < x\ < 0, this additional term will be significantly 
larger than zero only in the immediate vicinity of the trailing edge, provided that ek c C is large. 
Therefore, to include the effect of the incident pressure, the transfer function to be in Eq. (7a) is 
the two-component function g{x\, k Cl ^ 2 ) in [16], with k -2 = 0: 

g{xi,k c ,0) = e ek ^ - 1 + (1 + i)E*[- Xl {k c +n(l + M))} (7d) 
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where E* is the same complex combination of Fresnel integrals as in Eq. (3c). Amiet [5] was able to 
avoid the direct use of the parameter e because of the manner in which the transfer function in Eq. 
(7d) was used in his analysis. Amiet used the transfer function to define an unsteady lift response 
function that involved the chordwise integration of the transfer function with other terms. This 
integration yields an expression that, after clever manipulation of limits, does not contain e but still 
provides an additional term to the lift response function that accounts for the incident pressure. 
In the present case, the transfer function in Eq. (7d) must be explicitly used and therefore a value 
for e must be specified; a value is specified in the following subsection. 

To illustrate the effect of including this incident pressure term in the surface pressure formu- 
lation, Fig. 8 shows the same single frequency surface pressure fig. 3, with and without the the 
incident pressure term, with e = 1. Clearly, the incident pressure term has a significant effect only 
near the trailing edge, as expected. However, for a given frequency, differing values of e will affect 
differing amounts of upstream chordlength. Note that the addition of this term causes the pressure 
jump to vanish at the trailing edge for all time, i.e., the Kutta condition is satisfied. Note also the 
increase in spatial oscillation near the trailing edge when this incident pressure term is included. 

Eqs. (7a)-(7d) represent the the complex-valued broadband surface pressure formulation used 
for the present TE noise predictions. The final representation for the unsteady broadband pressure 
on the airfoil surface is then given by the real part of — A P in Eq. (7a). Using symmetry arguments 
and algebraic manipulation, the indicial bounds for the surface pressure spectral representation are 
altered so that the domain includes only positive wave numbers. As input to Eq. (1), the resulting 
real-valued surface pressure can be written 

N 

p(yi,r) = -47 r ^2 A n{ B n cos [k c ,n{yi - U C T ) + <f) n \ + D n sin[A: C! „(yi - U c t) + <p n \ } (8a) 

n= 1 
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where 


An = 


— h(u n ) S qq (u n , 0) A k c 

7T 


Bn = e^ n y. _ l+c ^ n) + s ^ n) 


D n = C(( n ) - S(( n ) 


£n — Vl[kc,n T /^n(l + M)\ 


(8b) 

(8c) 

(8d) 

(8e) 


and C(£ n ) and S{^ n ) are the Fresnel cosine and sine integrals in Eq. (3c). Specific evaluations for 
the correlation lengths I 2 and surface pressure correlations S« will be discussed in the following 
subsection. Furthermore, the zeroth term in Eq. (8a) is omitted because of the singularity in the 
formulation for the correlation length at n = 0. 


4.3. TIME DOMAIN PREDICTIONS 

The lower and upper frequency bounds for both calculations are 25 Hz and 10 kHz, respectively. 
Therefore, / = 25 Hz also serves as the fundamental frequency and the numerical bandwidth A/. 
Each calculation is performed for one period of the lowest frequency, T = 0.04 sec. The numerical 
solution is sampled at the Nyquist frequency, i.e. , At = T/2N. The calculation is performed on a 
500 x 100 surface grid with grid-point clustering near the trailing edge, as shown in Fig. 9. This 
trailing edge clustering allows for better resolution to account for the effect of the incident pressure 
term in Eq. (7d), as previously shown in Fig. 8. As in the constant frequency case, because the 
surface pressure in Eqs. (8a)-(8e) is cast in only one spatial variable y 1 , and the observer location is 
symmetric relative to the airfoil span, the acoustic predictions are found to be relatively insensitive 
to the discretization in 3 / 2 ! the primary concern for grid resolution is in the streamwise direction. 
With 500 points in the streamwise direction and clustering near the trailing edge, a sufficient 
resolution of at least 10 points per wavelength was obtained for the entire length of the chord. This 
conclusion was reached by inspecting surface pressure profiles for the highest frequency of 10 kHz. 

The coordinate system for the calculation is such that the xy-axis coincides with the center span 
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line, and the microphone position is in the plane X 2 = 0. The experimental microphone position for 
which comparisons are made is at a distance of 1.22 m from the model, at an angle of 90 degrees 
relative to the chord and directly above the trailing edge. Thus, the measured observer position 
for the prediction is, x = [0, 0, 1.22 ] T in meters. 

The baseline prediction case under consideration is for a tunnel speed of U = 69.5 m/s. This 
flow condition, the above observer location, and airfoil geometry are incorporated into an acoustic 
prediction using Eq. (1) with the surface pressure defined by Eqs. (8a)-(8e). Following the exam- 
ple of Schlinker and Amiet [16], a first-cut prediction is performed by using flat plate theory to 
evaluate the surface pressure correlations s m ( co,0) and spanwise correlation lengths in Eq. 

(8b). Schlinker and Amiet used empirical formulations for these quantities that they derived from 
previous analysis and boundary layer measurements of Corcos [13] and Willmarth and Roos [14]. 
The surface pressure correlations are approximated by 


S qq (w, 0) ~ q o 


<5* 2 x 1(T 5 

U 1 +w + 0.217 Co 2 + 0.00562 w 4 


(9a) 


where qo = p^U 2 / 2. 8* is the trailing edge displacement thickness, and Co = co8* jU . The displace- 
ment thickness is also taken from a flat plate approximation for turbulent boundary layer thickness 
8 on a flat plate, based on the chord Reynolds number Rec , he., 


8 0.37 5* 

O 

c (Rec)s C 


(9b) 


For the experiment of Brooks and Hodgson [17], Schlinker and Amiet [16] used Eq. (9b) to 
compute the boundary layer thickness 8 and accounted for boundary layer tripping by taking the 
15 percent chord station as the initial point of the calculation. Surface curvature was also accounted 
for in the downstream distance used in the calculation. The ratio 8 jC used by Schlinker and Amiet 
for this experiment was reported as 0.0166 for U = 69.5 nr/s and 0.0187 for U = 38.6 nr/s. The 
displacement thickness was then taken as 1/8 of the boundary layer thickness. The expression that 
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Schlinker and Amiet [16] suggest for the spanwise correlation length is 




2.H7 C 

00 


(9c) 


Fig. 10 shows the far field signal p'(x,t) that is predicted by Formulation IB at the experimental 
microphone location, for a tunnel speed of 69.5 m/s. The surface pressure is modeled with Eqs. 
(8a)-(8e) and (9a)-(9c). The time signal p'(x,t) is Fourier analyzed to determine a discrete set of 
spectral amplitudes {P n }n=i- The far field sound pressure level (SPL) spectrum is calculated by 


SPL(f n ) 



n = 1, 2, . . . , N 


( 10 ) 


where the reference pressure is P re f = 20/rPa. The SPLs are converted to a 1.0 Hz bandwidth by 
reducing the values in Eq. (10) by 10 log(A/). 

The resulting narrowband SPLs are compared with the prediction of Schlinker and Amiet [16] in 
Fig. 11. Also on this plot are the narrowband SPLs that were experimentally measured by Brooks 
and Hodgson [17]. The predicted results of Schlinker and Amiet and the measurements of Brooks 
and Hodgson shown in Fig. 11 were obtained by digitizing the appropriate plots of Schlinker and 
Amiet (Figure 34 of Reference 16). Various values of the parameter e in Eq. (7d) were assessed in 
this comparison stage of the research. The value e = 1.5 was chosen for its agreement with the flat 
plate correlation results of Schlinker and Amiet [16]. This value of e is held fixed for the remaining 
calculations. 

Clearly, Fig. 11 shows that significant error exists between the predictions and the measurements 
when flat plate formulations are used for the required surface pressure correlations. Fig. 12 sheds 
light on this error by comparing the flat plate formula in Eq. (9a) with the measured surface pressure 
correlations of Yu and Joshi [15]. The notation s qq denotes that the surface pressure correlations 
are normalized by 5*/U. The measured data in Fig. 12 were obtained by digitizing the “average” 
correlation plot of Schlinker and Amiet (Figure 35(a) of Reference 16). The normalized flat plate 
surface pressure correlations are significantly lower than the measured data, by as much as 7 dB. 
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The significant error in the flat plate approximation is only partially due the the lack of pressure 
gradient. The more significant error made in the approximation in Eq. (9a) is the lack of a trailing 
edge in its derivation. This empirical formulation is based on experimental measurement and 
analysis in which the flat plate is assumed to be infinite. Clearly, surface pressure correlations that 
are based on flat plate theory are inappropriate for predicting TE noise in this case. 

The experimental surface pressure correlations in Fig. 12 will now be used in the surface pressure 
formulation to predict the TE noise associated with the two tunnel speeds of interest and compared 
with experimental measurements. The modified formula for the surface pressure correlations is 

Sqq{v, 0 ) ~ Sqq(u,0) ( 11 ) 

where Sqq( tu. 0) denotes the normalized measured data in Fig. 12. The tabulated data for this mea- 
sured surface pressure correlation function is stored in a file that can be accessed and interpolated 
to obtain Sqq( cn, 0) for any frequency. The alteration of the surface pressure correlation function 
by experimental data now brings into question the evaluation of the spanwise correlation length. 
However, it was concluded by Brooks and Hodgson [17] that the function ^(w) for a flat plate and 
a thin airfoil are identical under suitable normalization. Therefore, the use of Eq. (9c) for ^(w) 
will be retained for the remaining calculations. 

The predicted and measured far field SPLs for the two tunnel speeds are shown in Fig. 13. 
The experimental data in Fig. 13 were obtained by digitizing the measurements plotted in Figure 
34 of Reference 16. The agreement with the measured data is significantly improved when the 
calculation includes surface pressure correlations that account for the trailing edge of an airfoil. In 
fact, it was concluded by Schlinker and Amiet [16] that airfoil surface pressure correlations were 
absolutely necessary for realistic TE noise predictions. 


CONCLUDING REMARKS 


The prediction of broadband trailing edge noise from rotating machinery and airframes is cur- 
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rently the subject of intense research in aeroacoustics. The physics of broadband noise generation 
are well understood as the result of the pioneering research of Howe [1, 24, 25], Amiet and coworkers 
[4, 5, 12, 16, 22], and Brooks and coworkers [17, 26]. The previous work of these and many other 
aeroacousticians has clearly demonstrated that any successful broadband loading noise prediction 
requires an understanding of two physical processes: the character of the time-dependent surface 
pressure that provides the acoustic source, and the manner in which that source gives rise to an 
acoustic signal. 

Obtaining the fluctuating surface pressure distribution analytically, numerically, or experi- 
mentally is itself a difficult problem. For this reason, past researchers have most often resorted 
to modeling the surface pressure, using guidance from experiments. However, as computational 
technology continues to advance, high resolution surface pressure fluctuations could soon become 
available from simulations of realistic problems. Therefore, the improvement of the acoustic radia- 
tion model becomes an important research topic. 

In the present work, a a loading- noise solution of the Ffowcs Williams-Hawkings equation is 
applied to the prediction of broadband trailing edge noise on a NACA 0012 airfoil in a low Mach 
number flow. This new solution, originally proposed in [7], is to date the simplest analytical result 
for the prediction of loading noise and is suitable for statistical analysis of broadband noise for a 
surface in general motion. The time domain predictions are found to be in excellent agreement 
with the frequency domain predictions of Schlinker and Amiet [16] as well as with the experimental 
measurements of Brooks and Hodgson [17]. These results are, to the authors’ knowledge, the first 
successful broadband trailing edge noise predictions in the time domain. 
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APPENDIX: LIST OF SYMBOLS 


b 

C 

co 

/ 

/ 

E* 

9 

k c 

h 

M 
M r 
M v 
A P 
P 
P' 

dp/ ds 

9o 

r 

s qq 

T 

t 

U 

V 

x 

y 


airfoil semi-span (m) 
airfoil chord (m) 
ambient sound speed (m/sec) 
frequency (Hz) 

geometry function for airfoil surface (Fig. 1) 
combination of Fresnel integrals (Eq. (3c)) 
surface pressure transfer function 
<x/U c . convective wave number (rad/m) 
spanwise correlation length (m) 

P/co, Mach number vector 

M ■ r/r. Mach number in radiation direction 

M ■ Mach number in direction of u 

unsteady airfoil surface pressure jump (Pa) 

unsteady airfoil surface pressure (Pa) 

sound pressure radiated to observer (Pa) 

surface pressure gradient in the direction of V (Pa/m) 

PqU 2 /2, dynamic head (Pa) 

x — y , sound radiation vector (m) 

surface pressure correlation function (Pa 2 -sec/rad 2 ) 

1 // , acoustic period (sec) 
observer time (sec) 
uniform freestream speed (m/sec) 
airfoil velocity vector (m/sec) 
pi, X 2 , x^Y i observer position (m) 

[yi,?/ 2 , 0 ] T , surface source position (m) 
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p = Vl - M 2 

6 = boundary layer thickness 

5* = displacement thickness 

A = Co//, acoustic wave-length (m) 

p = Moj//3 2 U (rad/m) 

ip = directivity angle (Fig. 5) (rad) 

v = unit inward facing normal on surface edge (Fig. 1) 

6 = angle between surface normal and r (Fig. 1) (rad) 

Pq = ambient density (kg/m 3 ) 
r = source time (sec) 

<IVp = power spectral density of surface pressure (Pa 2 -m 2 /rad 2 ) 
ip = random phase variable (rad) 
oj = 27 r/, circular frequency (rad/sec) 

Subscripts 

1,2,3 = Cartesian coordinate directions 
ret = evaluated at retarded time t — r /cq 
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Figure 2 Schematic for the constant-frequency trailing edge noise problem in Section 3. 
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Figure 4 Schematic for directivity calculation, observer on circular path in plane X 2 = 0. 
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Figure 5 Directivity for a constant frequency source of 2.5 kHz; observer distance 2 m; spacing 
between concentric circles on grid represents 0.25 Pa. 
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Figure 6 Velocity scaling properties as determined by Formulation IB and the surface pressure 
in Eqs. (3a)-(3d). — Q — , prediction; , slope = 5. 
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Figure 7 Schematic for trailing edge noise experiment of Brooks and Hodgson [17] 
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Figure 8 Effect of incident pressure term on initial surface pressure profile using Eqs. (3a)-(3d), 
and modified with Eq. (7d). , induced pressure; , sum of induced and incident pressure. 
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Figure 11 Predicted and measured far field noise spectra; predictions performed with surface 
pressure correlations from flat plate theory; U = 69.5 m/s. O • frequency domain prediction 
from [16]; , time domain domain prediction; , experimental measurement from [17]. 
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Figure 12 Normalized surface pressure correlations. 
— O — j experimental measurements from [15]. 


, from flat plate theory (Eq. (9a)); 
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Figure 13 Predicted and measured far field noise spectra; predictions derived from time domain 
formulation (Eq. 6) and measured surface pressure correlations [15]. A , prediction for U = 38.6 
m/s; O 5 prediction for U = 69.5 m/s; , experimental SPLs from [17]. 
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